SOLUTION OF TWIST-3 EVOLUTION EQUATION IN 
DOUBLE LOGARITHMIC APPROXIMATION IN QCD 



A.G.Shuvaev 

Petersburg Nuclear Physics Institute, 
Gatchina, St. Petersburg 188300 Russia 



Abstract 

The solution of DGLAP evolution equation for the twist-3 gluon 
operators is obtained in the Double Logarithmic Approximation of 
QCD perturbation theory. The method used for the solution is similar 
to the reggeon field theory. The asymptotics of the twist-3 parton 
correlation function for the small Bjorken variables xb is found. 

1. The growth of the structure function f(xB,Q 2 ) in the region of small 
Bjorken variables xb makes it necessary to take into account the high twist 
contribution. According to the operator product expansion the local oper- 
ator of spin J contributes to the number J Mellin moment of the structure 
function f(xB,Q 2 )- The Q 2 -dependence of the operators is perturbatively 
determined by evolution equation collecting in the leading order (LLA) the 
terms of the form (as lnQ 2 //i 2 ) n . The asymptotic behavior for xb — > is 
governed by the rightmost singularity of the anomalous dimension in the 
variable J continued to the complex plane. However the calculation of the 
anomalous dimension becomes more complicated for twists N > 3 because 
of the large number of the local operators. 

There is another approach proposed in Ref.fl] and based on the relation 
of small-xs behavior of the structure function and BFKL equation summing 
up the powers (agin 1/xb) 71 - It enables to find twist-2 anomalous dimension 
near the singularity position J — > 1 for the leading and next to leadings 
orders in terms of InQ 2 . The situation is also more involved for the higher 
twists N > 3 since one has to solve the equation for N reggeized gluons [2] 
and then to extract from it the anomalous dimension. 

2. Here we consider the evolution equation for the twist-3 quasipartonic 
operators in double logarithmic approximation (DLA), which collects the 
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powers of the product as lnQ 2 //! 2 In 1/x. Quasipartonic operator form the 
closed set of high twist operators allowing for the interpretation in terms of 
the parton model [5]. They are responsible for the small xb asymptotics of 
the structure function [3, 4]. The matrix elements of quasipartonic operators 
depend only on the fraction x,i of the partons momenta along the hadron 
momentum p (p 2 ~ 0). The pure gluon channel will be studied below as 
dominating in the small Xb region. We shall take the quasipartonic operators, 
which in the axial gauge n^A^ = with a light-like vector n dual to the 
hadron momentum p, have a general form 

= r M tl^^((<9) mi ^)((<9)"«^)((<9)™-^), 

where m 8 are positive integers, T is a color and Lorentz tensor, the particular 
form of which will not be important in what follows, d = n^d^. The matrix 
element over a hadron state can be expressed as 

/n abc a mi a m 2 a m 3 

x Ly Ai/ii,A 2 /i2,A 3 ^3^1 ^2 P3 > 
/nabc _ -pafec ^.Ai A 2 A3 

Ly Ai/ii,A 2 /i2,A 3 ^3 ~~ 1 /ji/i' 1 ,/j2M 2 >A»3A»3 Mi M 2 M 3 ' 

where is the gluon polarization vector. The parton correlation function 
N^x 2 x 3 (xi, xii X3) has a meaning of a hadron wavefunction integrated over 
partons transverse momenta, the greatest transverse momentum being of the 
order Q 2 . 

The evolution equation generalizing twist- 2 DGLAP equation is derived 
in paper [5] . It has the form of iV-particle one dimensional Schrodinger-type 
equation with pairwise interaction between gluons, 

r v) ^/ (i) 



i<j 

ai,...,a'a'...,a N 



X ^X i \ J ,\' i \'i X ^ X j'^i^j)^X 1 ,...,X' i ,X' j ,...,X N (Q ' X l' • • • ' A) • • • ) Pji ■ ■ ■ 

This equation sums up in the leading \nQ 2 order the ladder-type diagrams. 
For the twist iV case they comprise the local operator vertex and N gluons in 
t-channel interacting through all possible s-channel gluons rungs. The inte- 
grals in each ladder cell are ordered in LLA such that transverse momentum 
in the above cell plays the role of ultraviolet cut-off for the below one. The 
Q 2 value, being the greatest momentum in the upper loop attached to the 
local operator vertex, is the ultraviolet cut-off for the whole diagram. The 
evolution equation is obtained by taking the derivative of the diagrams with 

respect to logQ 2 . The kernel (xj, Xj\ Pi, fy) is determined by the log- 

arithmic part of the one-loop integral over the transverse parton momentum 



Generally the longitudinal momenta Xj are not ordered in LLA, but they 
have to be ordered in DLA to provide a large logarithm for each ladder cell. 
The Xi variables increase from the smallest values at the local operator vertex 
to the order of unity ones in the lower part of a diagram. In such a kinematics 
the logarithmic divergencies that occur in every loop when fa — > is cut from 
below by the longitudinal momentum in the upper cell. Thus DLA implies 
the loop integrals in the evolution equation (1) to be limited by the condition 

Pi, Pj Xi,Xj, (2) 

which means that the momenta below the s-channel rung (x) and above it {fa) 
are of different order of magnitude. The most singular contribution comes in 
Eq.(l) from the region where both Pi and Pj tend to zero. The momentum 
conservation allows it only if 

Xi + Xj < Pi,Pj, (3) 

that with logarithmic accuracy. Hereafter it is convenient to 

assume the momenta directed upward to be positive, directed downward to 
be negative. 

In the logarithmic domain the kernel of the evolution equation in DLA 
can be easily obtained by keeping the most singular in ft terms in the gluon- 
gluon kernel, 

<: c xU4 x ^ x ^^^ = ^M^r H f b9d x ^ ~ X2) j[p 2 (4) 

Here \, A- are the two dimensional transverse helicity indices, f abc are the 
structure constants of the SU(N C ) group. The momenta x±,x 2 , are equal 
in DLA. They are positive but have the opposite directions, one of them is 
incoming from below the other is outgoing. The low-scale momenta Pi j2 are 
not supposed to be equal in DLA since the momentum transfer from below 
Xi — x 2 is small only compared to the large momenta x i;2 but is of the same 
order as the low-scale ones. 

All Xi, Pi momenta are supposed to be positive in the DLA kernel (4), 
the sign is specified with an additional index a = {+,—}. Thus each gluon 
in the structure function is characterized by the color index a, helicity A, 
momentum value P and momentum direction a. The interaction occurs only 
between the gluons with the opposite a. 

There are two possible color structures for twist-3 operators - the odderon- 
like d abc and gluon-like f abc . Both of them go through the equation resulting 
into N c /2 factor. This simplifies the color structure of the correlation func- 
tion, = ^ c F Al , A2iA3 or = f abc F XlMM , and the action of the 

kernel can be written as 

Hi2Fx u<Tl ,x 2t<r2 ,x 3 , <T3 (xi,X2,X3) = -a5 ai - a2 5 Xl ,\ 2 5\' 1 ,\> 2 x 1 5(x 1 - x 2 ) (5) 



f xi dpi f X2 dp2 

X Jo ~p\ Jo ~fc~[ Fx >uX 2 ,<r2,M,*3(f3l,f32,X 3 ) ± F x > 2!a2!X > vai! x 3t(73 (f32,Pl,X3,)}, 

- - AT aS 
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and similarly for H23, H13. Two terms in the r.h.s. (5) represent the sum 
of the s- and -u-channel diagrams (only even momentum operators survive in 
this sum for the twist- 2 case). 

3. Instead of direct solving of evolution equation we adopt here an other ap- 
proach similar to reggeon calculus and more suitable to find the asymptotics 
of the structure function. To this end we rewrite the formal solution of the 
evolution equation with a given initial condition F , 

F(Q 2 ) = e Hl °z Q2 ^ 2 F , 

through the Mellin transform as 
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where the integral runs along the imaginary axis to the right from all singu- 
larities. The equation 

F(u) = F + ^-H F(v) (6) 

can be treated as the Bethe-Solpeter equation in the theory described by 
effective action 

S = f dx<$> Xa (x)$ Xa (x) 

$>l z(J3 {x)6{x-(3 2 )<$> XAaA {(5 2 ). 
The solution to eq.(6) is given by the convolution with the Green function 

F x >a>(is, fr) = J dxitp^x^G^xy^Xirfi) (7) 
calculated in the effective theory, 
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1=1 

The initial hadron wavefunction tp\ iai ( x i) is to be taken at the low Q 2 scale. 
It can't be find perturbatively, but its precise form do not influence large Q 2 
asymptotics. For defmiteness we shall consider below the moments, that is 
we take 

<p(xi) = flx? (8) 

i=l 

with the integers rii > 0. 



Figure 1 



We shall find the Green function by iterating the Bethe-Solpeter equation. 
We start from the two gluon ladder, or "reggeon", which will be a main 
building block in the further proceeding. It is schematically shown in Fig.l, 
where the solid lines denote gluon. Iterations of the two-particle kernel results 
into the matrix 

1 11 

9x 1 x 2 ,x' 1 x' 2 (xi,x 2 ;pi,p 2 ) = - 5 x 1 \ 2 5x> 1 x>——gi2(v;xi,{{3 1 ,f3 2 }) (9) 

x x\ 5(xi — x 2 ). 

Here /3 2 } = max {Pi, P 2 } and matrix g ik acts on the sign variables at, a k 

as 

g ik (v;x,P) = g Cia ^ a y k {v;x,P) = ( I±Ptk A ik ^ g{y;x,p), 



k 



where operator P ik permutes the indices Oi, a k , matrix A ik permits the in- 
teraction only between the partons with opposite momentum signs, 



J _/ 



and 

4. The expression (9) leads to the usual structure function in the twist-2 
case. Indeed, the general form of twist-2 spin J gluon operator structure 
function) convoluted with the gauge fixing vector n is 

n ltl ---n (iJ (ilt ... tliJ = {id)A v {id) J ~ 1 A v . 

It results into the vertex 

5x 1 x 2 prp 2 n2 m+P2) 

with mi — 1, m 2 = J — 1. There is no momentum transfer through the 
operator, this is the reason for the momentum delta-function. This vertex 
should be integrated with the ladder function (where a/2 is replaced with a 



for the color singlet). The integration has to be done with account of both 
signs of the /?i ; 2 variables, which implies the sum over (H — ) and ( — h) initial 
sign states. For the (H — ) final state, that is for the positive x\, we get 



M 2 (J) = -[l + (-iy]5 XlX2 x{5(x 1 -x 2 ) J 



dj a 



2m v(J — 1) — a 



and the same for the ( — h) state. Thus we have reproduced the double 
logarithmic twist-2 anomalous dimension 



J -I 



together with the selection rule allowing only for the even J values. 
5. We consider the diagrams for the Green function occurring in the eq.(7). 
The general sum of the three-gluon ladder diagrams can be equivalently pre- 
sented as a sum of two-gluon ladders ("reggeons") developing between each 
gluon pairs accompanied with a third single gluon as is shown in Fig. 2 - 4. 
By employing this representation all diagrams can be summed up in a closed 
form. The Green function reeds 

G tf = £ Guixi^WPv^i'M (ii) 

where the sum is taken over independent permutations of the incoming and 
outgoing particles while the functions Gdj stand for the diagrams with a fixed 
order of the external lines. The symbols d and / label the Green functions for 
the d abc and f abc color structures. The formula (11) implies the simultaneous 
permutations of all quantum numbers, that is momenta, helicities and colors, 
which means symmetrization with respect {xi, <Jj, A«} (or {flj, a'j, Xj}) pair for 
d abc tensor and antisymmetrization for f abc tensor. 



xi,ai, Ai 




X2,&2, A 2 



$5, 03, ^3 



X3,cr 3 , A 3 



Figure 2 

The effective diagrams constructed from the two-gluon ladder and gluon 
line turns out to be rather simple to calculate the Green function by the 
direct summation. The result is presented by the sum of three contributions. 
Two of them are degenerate in the sense that they comes from the finite 
number of iterations. The first one includes the "reggeon" only once. There 
are 3 diagrams of this type for 3 various gluon pairs combined into the ladder. 



One of them for the pair 12 is shown in Fig. 2, 

1 1 

Gdj (x u A)(i2) = 2 5 a 1 a 2 5a 3 a^a;a^ j-^ £12 (f| ^i, {A* &}) 
x xiS(xi - x 2 )5(x 3 - (3 3 ) 8 a3a > 3 . 

The others can be obtained by permutations of the indices (123) — > (231) 
and (123) -> (132). 



A,ai,A' /3 2 , cr 2 , A' 2 



Xi,(Ti, Ai 




( 12 ) R rr> V 



(23) 



x 2 ,ct 2 , A 2 x \ x 3 ,(7 3 , A 3 
Figure 3 

The second contribution arises from the diagram with the two "reggeons" . 
Fig. 3 presents the diagram where the gluon pair 12 switches to the pair 23, 

Gaj (xi, A)(23)(i2) = ^ S X2 x 3 5 Xl x' 3 5 K y 2 p g 23 {v; x u {pi, f3 2 }) (12) 

x gi2(v,x 2 , {x!,/3 3 }) x 2 S(x 2 - x 3 ). 

The other 5 terms result in this case from (12) after independent permuta- 
tions (123) — > (231) and (123) — > (132) of the upper and lower (in the sense 
of Fig. 3) indices but excluding the equal ones. In other words the sum is 
taken over various ways to combine the incoming and outgoing gluons into 
different two-particle ladders. 

The contributions starting with the three "reggeons" develop a regular 
series which can be written as 

X 1 

Gaj (Xi,A)(12)(12) = ^ ^A^AsA^A'^ p ± p 2 p 3 g!2( y ; X U X3) (13) 



x J ^W d j(x 3 ,v;{(3 3 ,(3})g 12 (v;P,{Pi,fc}), 
W dJ (is;x,P) = \g(y;x,P) (14) 



+ 



\f^-g(u;x,(3')g(u;(3',(3) + 



Fig. 4 shows the diagram corresponding either to the first term in eq.(14) or 
to all terms in (14) if the middle "reggeon" is replaced with the full function 




Figure 4 

There are 8 other terms besides (13), which can be obtained from it by 
independent permutations (123) — > (231) and (123) — > (132) of the upper 
and lower indices in the Fig. 4. 

The Mellin transform (10) turns the convolutions over f3' into the usual 
products, making the series (14) to be equivalent to the purely matrix prob- 
lem, 

^ = 5^ + ^ , + - = ff 2?b- (15 » 

where the "reggeons" "dissociate" into two-particle interaction, 

H = H itk = -S Xi \ k Sx'.x' k — A ik , 

and the problem is reduced to the inversion of the finite matrix. 

The above iterations exhibit that at each step two momenta with oppo- 
site directions have the same value much larger than the value of the third 
momentum (#2,3 %i in Fig.3 and 3> X3 in Fig. 4). This property ex- 
preses the longitudinal momentum conservation within the DLA accuracy - 
the sum of all momenta is small compared to their natural scale. 
6. The Green function is convoluted over variables $ and helicities indices 
A^ with the operator vertex given by the expression 

(16) 

where the longitudinal (^-function corresponds to forward kinematics without 
momentum transfer. Taking then the moments with respect variables x« 
(8) we get the Green function in the moments representation, G^°j(mj, n^). 
Note, that the integrals over Xi, Pi imply the positive as well as the negative 
values of the momenta. The negative values are described through the sign 
variables <7j = ±, for example, the configuration where (5\ < 0, {3 2, 3 > is 
associated with the state (—,+,+) and similarly for x { . The integral over 



all sign configurations is recovered by the sum over all initial a[ and final 
Oi values. As a result the Green function written in terms of the moments 
takes into account both signs of $ and Xi and does not contain the auxiliary 
variables a iy 

{i'},{i} 1 2 3 

(17) 

The sum here means the independent symmetrization with respect the pairs 
{mj/,Aj/} and {rit, Aj} for d abc and antisymmetrization for f abc structures. 
The helicities A- should be convoluted with the tensor Fy. , w specifying the 
operator vertex (16). 

Separating the common factors the functions G d j(rrii, ni) takes the form 



Gufam) = [l + (-l) m+ 1((-in±(-l) 



m-2 mi 



X 7:G d f(mi,ni), (18) 

m + n + 2 ' A ' ; ' v ; 

m = m l + m 2 + m 3 , n = n x + n 2 + n 3 , 

where + and — stand for d and / structures, respectively, and the function 
Gd,f is expressed in terms of the three contributions considered above, 

The first contribution yields 

G S( m "^ = l2^i' (19) 
j = m + n 3 , 
while the second one takes the form: 

G d n \m t ,n t ) = - , * ^— (20) 

rf v " 11 8 2i/(m - 1) - a 2vj - a K J 

Gf ] = 0. (21) 

(G^ vanishes after antisymetrization over end points). The third contribu- 
tion with the matrix (15) inverted reads 

c (ni) 3 a a a 

d B2u(m-l)-a2uj -a 4i/(m-l)-3a' 1 ' 

G (iii) 3 g « g / 23 n 

f ~ 8 2z/(m- 1) -a 2vj -a Av(m - 1) -a' 

7. The asymptotic behavior of the structure function for Bjorken variable 
x B — > is determined by the rightmost singularity in the variable J, which 



has the meaning of the local operator spin continued to the complex plane. 
The spin of quasipartonic operator J — m, therefore one needs to continue the 
function Gdj(rrii, rii) to m — > 1 formally keeping the other variables rrii, rii to 
be fixed. Because of the signature-like factors (— l) m+ra the terms with even 
or odd m are to be treated separately. Note that in a general LLA case this 
continuation is non trivial since the mixing matrix describing the evolution 
has a rank depending on J [8]. The explicit form of the DLA solutions (19) 
- (23) makes the continuation much more simple and straightforward. The 
obtained results show that the anomalous dimension for d abc color structure 
is _ 

^ = i7^T- < 24 > 

The main singularity for f abc structure is actually given by the pole of 2 
gluon state, 

7/M = lj^~ r (25) 

The contribution of " developed" 3 gluon ladder is more weak in this channel, 
7j(J) = 7/(J)/2. The contributions of other singularities are strictly speak- 
ing beyond the DLA accuracy since they produce the extra positive powers 
of Xb- 

The DLA anomalous dimension (24) is smaller compared to those, which 
can be derived from the direct solution of BFKL equation obtained in refs.[4, 
6, 7, 9]. The possible reason for this is that the known BFKL solutions 
founded for Odderon do not really correspond to quasipartonic operators of 
twist 3. In this case the DLA result (24) could indicate to an existence of 
another solution of quasipartonic type. 
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